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Let U n be an n x n Haar unitary matrix. In this paper, the asymptotic 
normality and independence of Tr U n , Tr U%, . . . , Tr are shown by us- 
ing elementary methods. More generally, it is shown that the renormalized 
truncated Haar unitaries converge to a Gaussian random matrix in distri- 
bution. 
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o 1 Introduction 



Entries of a random matrix are random variables but a random matrix is equivalently 
considered as a probability measure on the set of matrices. A simple example of random 
matrix has independent identically distributed entries. In this paper random unitary 
matrices are studied whose entries must be correlated. 

A unitary matrix U = (Uij) is a matrix with complex entries and UU* = U*U = I. 
In terms the entries these relations mean that 

n n 

Yl i^i 2 = 2 i^i 2 = ^ ioT a11 1 - ^' - n ' w 

j=l i=l 
n 

U a U jh for all 1 < i, j < n, i ^ j. (2) 

i=i 
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The set Uin) of n x n unitary matrices forms a compact topological group with 
respect to the matrix multiplication and the usual topology, therefore there exists a 
unique (up to the scalar multiplication) translation invariant measure on U(n), the 
so-called Haar measure. We will consider a random variable U n which maps from a 
probability space to U(n), and take its values uniformly from U(n), i.e. if H C U(n), 
then 

Prob (U n eH)= 1 (H), 

where 7 is the normalized Haar measure on U(n). We call this random variable a Haar 
unitary random variable, or shortly Haar unitary. 

Relations (1) and (2) show that the column vectors of a unitary are pairwise oth- 
ogonal unit vectors and the distribution of each column has unitarily invariant joint 
distribution in case of a Haar distributed unitary. This fact allows a simple construction 
of a Haar unitary from a Gaussian matrix with i.i.d. entries. In Sect. 2 this construc- 
tion is written up in details. The construction allows to determine the distribution of 
the matrix elements. Due to the permutation invariance, the elements are identically 
distributed. We observe that large powers of the eigenvalues are independent and uni- 
formly distributed. The main aim of the paper is to study asymptotic questions when 
the matrix size is going to infinity. In this limit the matrix elements are going to be 
Gaussian (after a renormalization) and more generally, the truncated matrix converges 
to a Gaussian matrix. This is the content of Sect. 3. In the rest of the paper we show 
that the trace of the powers is going to Gaussian in the limit, moreover the traces 
of different powers are asymptotically independent. Actually, this has been shown by 
Diaconis and Shahshahani [3] , they determined the Fourier transform of the limit dis- 
tribution of the eigenvalues. In their proof the characters of the symmetric and unitary 
groups, and different bases of the symmetric polynomials play the main role. Here we 
get the same result in a more elementary way. The method of moments is used, and we 
examine the order of magnitude of the different summands in the trace. This method 
could be familiar from Arnold pQ when he studied Wigner matrices, or from Bai and 
Yin [2J for sample covariance matrices. 

Random unitary matrices may be applied to several physical phenomena, such as 
chaotic scattering or statistical properties of periodically driven quantum systems [Sj. 
The large deviation theorem for the emprical eigenvalue density of Haar unitaries (and 
of some other unitaries) was established in [7j. Independent Haar unitaries provide also 
a simple example of asymptotical freeness [SJ El • 

2 Haar unitary matrices 

Let £ be a complex- valued random variable. If Re£ and Im£ are independent and 
normally distributed according to N(0, 1/2), then we call £ a standard complex normal 
variable. The terminology is justified by the properties E(£) = and = 1. 
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Lemma 2.1 Assume that R > and R 2 has exponential distribution with parameter 
1, 9 is uniform on the interval [0, 2tt], and assume that R and 6 are independent. Then 
£ = Re ld is a standard complex normal random variable and 

E(ef)=hek\ (k,£eZ + ) 

Proof. Let X and Y be real-valued random variables and assume that X + iY is 
standard complex normal. For r > and < 9 < 2n set 

S rfi0 := {pe^:0<p<r,0<^<6 }, 



then 



P{X + iYeS rfio ) 



7Y 




2+t2 )dsdt 



' {(s,t):s+iteS rA) } 

i r e " r 2 

= - / # / pe~ p dp 
x Jo Jo 

= "U> (l-^ r2 ) =P(£eS r ,e ). 
This proves the first part which makes easy to compute the moments: 
E{^ k f ) = E (R k+i ) E ( e w ^) = 5 M E{R 2k ), 

so we need the moments of the exponentional distribution. We have by partial integra- 
tion 



f 

Jo 



x k e x 



k— 1 „— x 



e x dx = k 



' dx = - [x k e x ] ™ + k 

Jo 

poo r 

= fc(Jfe-l) / x k ~ 2 e~ x dx = . . . = k\ 

Jo Jo 

which completes the proof of the lemma. 



x 



k—l„—x 



dx 



e x dx = k\ 



□ 



Next we recall how to get a Haar unitary from a Gaussian matrix with independent 
entries by the Gram-Schmidt orthogonalization procedure on the column vectors. 

Suppose that we have a complex random matrix Z whose entries Z^ are mutually in- 
dependent standard complex normal random variables. We perform the Gram-Schmidt 
ortogonalization procedure on the column vectors Zi {i — 1,2, ... ,n), i.e. 

i-l 

Ui = ^ , (3) 



i-l 



where 



||(Xi,X 2 , . . . ,X n ) 



\ k=l 
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Lemma 2.2 The above column vectors U, constitute a unitary matrix U = (Uj) i=l . 
Moreover, for all V 6%) the distributions of U and VU are the same. 

Proof. First we prove, that for any V 6%) the matrices Z and VZ have the same 
distribution. The entries £y of VZ are standard complex normal. Indeed, 



i=i 



is normal. Furthermore 



£(&) = Ews«) = i), 



and 



Next we prove that the correlation between two entries is zero. (In the case of nor- 
mally distributed random variables this is equivalent to the independence.) 



mjU = E [[J2 VuZ y (£ VskZk 

n n n 
= E E skE(ZijZ kr ) = 5j r E VilVsl = 5j r 5 is . 



1=1 k=l 



1=1 



The ith column of VU is exactly VUi and we have 



i-l 



i-1 



VZt - ^{Zi, UfiVUt VZ, - J2( VZ i, VUiWUt 
\ 1 . = ^ = ^ (4) 



i-l 



z t -J2(z t ,u l }u l 



1=1 



VZi-Y.fyZ^VU^VUt 



1=1 



which is the Gram-Schmidt ortogonalization of the vectors VZ,. Since we showed above 
that Z and VZ are identically distributed, we conclude that U and VU are identically 
distributed as well. Since the left invariance characterizes the Haar measure on a com- 
pact group, the above constructed U is Haar distributed and its distribution is right 
invariant as well. □ 

The column vectors of a unitary matrix are pairwise orthogonal unit vectors. On 
the bases of this fact we can determine a Haar unitary in a slightly different way. The 
complex unit vectors form a compact space on which the unitary group acts transi- 
tively. Therefore, there exist a unique probability measure invariant under the action. 
Let us call this measure uniform. To determine a Haar unitary, we choose the first 
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column vector U\ uniformly from the space of n- vectors. U 2 should be taken from the 
n — 1 dimensional subspace orthogonal to U\ and choose it uniformly again. In general, 
if already Z/i, U 2 , ■ ■ ■ , Uj is chosen, we take Uj + i from the n — j dimensional subspace 
orthogonal to U±, U 2 , ■ ■ ■ , Uj, again uniformly. The column vectors constitute a unitary 
matrix and we check that its distribution is left invariant. Let V be a fixed unitary. 
We show that the vectors VUi, VU 2 , ■ ■ ■ , VU n are produced by the above described 
procedure. They are obviously pairwise othogonal unit vectors. VU\ is uniformly dis- 
tributed by the invariance property of the distribution of U\. Let V(l) be such a 
unitary that V(l)VC/i = VU X . Then V^V^VUi = U x and the choice of U 2 gives that 
V- 1 V(1)VU 2 ~ U 2 . It follows that V(1)VU 2 ~ VU 2 . Since V(l) was arbitrary VU 2 
is uniformly distributed in the subspace orthogonal to VU\. Similar argument works 
for VU3, . . . , VU n . The Gram-Schmidt orthogonalization of the columns of a Gaussian 
matrix gives a concrete realization of this procedure. 

The permutation matrices are in U(n), and by multiplying with an appropriate 
permutation matrix every row and column can be transformed to any other row or 
column, so the translation invariance of a Haar unitary U implies that all the entries 
have the same distribution. From the above construction of a Haar unitary one can 
deduce easily the distribution of the entries: 



\fnUij converges to a standard complex normal variable. The correlation coefficient 
between \Ua\ 2 and \Ujj\ 2 is l/(n — l) 2 if i ^ j (see p. 139 in 6J). 

In the next section we need the following technical lemma which tells us that the 
expectation of many product of the entries are vanishing. 

Lemma 2.3 Let i u . . . , i h , j u . . . j h e {1, ...,n} and k u . . . , k h , m 1 , . . . , m h be 

positive integers for some h G N. If 




n 



- 1 



71 



(see also p. 140 in 0). Since 





for some 



1 < u < n 



or 




for some 



1 < v < n. 



1 ) 



then 
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Proof. Suppose that t := J2 ir=u (K — m r ) ^ 0. The translation invariance of U implies 
that multiplying this matrix by V = Diag(l, . . . , 1, e l9 , 1, . . . , 1) e U(n) from the left 
we get 

E {iU^JJ^J ■ ■ ■ iU^UZj) = e U6 E ((U^UT^) . . . {U^ jh K)S) > 

for all 9 e R. □ 

Let C/ be a Haar distributed n x n unitary matrix with eigenvalues A , Ai, . . . , A n _i. 
The eigenvalues are random variables with values in T := {2 6 C : \z\ = 1}, their joint 
distribution is well-known: 



jfflll*-*l'=slll«"-«*f © 

l<j Kj 

with respect to dz dz\ . . . dz n _i, where dzi = d9i/ 2tt for Z; L = e ldi (see p. 135 in 6J). 



Theorem 2.4 For m > n the random variables A™, A™, . . . , A™_ 1 are independent and 
uniformly distributed on T. 

Proof. Since the Fourier transform determines the joint distribution measure of A™, A™, 
• • • , K1-1 uniquely, it suffices to show that 

f z k rz k r . . . zt~r n 1* - *,f <** = (6) 

j i<j 

if at least one kj G Z is different from (dz = dz dz\ . . . dz n _i and integration is over 
Let 

A(*b, zi, . . . , z„_i) := JJ(^ - = det[^ fc ] < i < n „ 10 < fc < n _ 1 . (7) 

(What we have here is the so-called Vandermonde determinant.) Then 
[ \ \%i — z j\ 2 — A(zb, 2Ti, . . . , z n _i)A(^ 1 , z x 1 , . . . , z n \) 

i<j 

and one can write (jHJ) as an n-times complex contour integral along the positively 
oriented T: 

Zq Z-y ... Z n _-y IXyZQ, Zi, . . . , 1 )A( 

= / ^ om • • • E (-i) ctW 4 (0) • • • ir 1} E (-ir (p v (0) ■ ■ ■ *£r l) ^ 

which is calculated according to the theorem of residue. We need to find the coefficient 
of Zq 1 z^ 1 . . . Zn-i, so we are looking for the permutations, for which kjm+ir(j) — p(j) = 
— 1 if < j < n — 1, so kjin = p(j) — 7r(j) — 1. Here \p(j) — n(J)\ < n — 1, and 
\kjm\ > m > n, if k 3 • ^ 0, so if at least one fcj e Z is different from 0, then there exists 
no solution. This proves the theorem. □ 
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3 Asymptotics of the trace 



Let U(n) = (U(n)ij) : Q — > U{n) be a Haar distributed unitary random matrix. In this 
section we are interested in the convergence of Tr U (n) as n — > oo. Since the correlation 
between the diagonal entries decreases with n, one expects on the basis of the central 
limit theorem, that the limit of the trace has complex normal distribution. We prove 
this by the method of moments. 

Theorem 3.1 Let U(n) be a sequence of n x n Haar unitary random matrices. Then 
Tr U(n) converges in distribution to a standard complex normal random variable as 
n — > oo. 

Proof. For the sake of simplicity we write U instead of U(n). First we study the 
asymptotics of the moments 

n , n , 

E((TrU) k (TFU) k )=E^J2 U ^ (E 17 ^) ) 

i=l j=l 

n n 

= EiJJi^ . . . Ui k i k Uj t ji ■ ■ ■ Uj k j k ), 

i 1> ...,i k =ljx,...,j k =l 

k e Z + . By Lemma \'2. 31 parts of the above sum are zero, we need to consider only those 
sets of indices {ii, . . . , i^} and {ji, . . . ,jk} which coincide (with multiplicities). Look at 
a summand E(\Ui 1 i 1 \ 2kl . . . \Ui r i r \ 2kr ), where YTi=i ^ = ^- ^F rom the Holder inequality 

EdU^r ■ ■ ■ \U irj f kr ) < II *{/ E{\U kjl \^) = n ( U H % 7 = O (n- k ) . 

1=1 1=1 ^ ' ' 

(8) 

The number of those sets of indices, where among the numbers ii, . . . , there are at 
least two equal is at most 

jfe!^7i fc - 1 = 0(n fc - 1 ). 

By (|HJ) the order of magnitude of these factors is 0(n~ k ), so this part of the sum tends 
to zero as n — > oo. 

Next we assume that ii, . . . , are different. Since by translation invariance any row 
or column can be replaced by any other, we have 

E(\U tlll \ 2 . . . \U ikik \ 2 ) = £(|f/ n r • • • |^| 2 )) =: Ml. (9) 

It is enough to determine this quantity and to count how many of these terms are in 
the trace. 



7 



The length of the row vectors of the unitary matrix is 1, hence 



J2--T l E(\U ill f l ...\U ikh \ 9 ) = l. (10) 

'•1=1 ifc=l 

We divide the sum into two parts: the number of terms with different indices is 
n!/(n — k)\, and again the translation invariance implies that each of them equals 
to M£, and we denote by e k the sum of the other terms. Therefore 



* = 1 -(^^* l (*, ,of,r * wo - 



and 

M n (1 ~ e n k )(n - k)\ 
k n\ 

Now we can count how many expectations of value M k are there in the sum (JHJ). We 
can fix the indices i\, . . . ,i k in n\/(n — k)\ ways, and we can permute them in k\ ways 
to get the indices ji, . . . , j k . The obtained equation 

lim E ((Tr^) fc (Trf4) fc ) = lim — ^l— fc! (1 " e£) ^ " = k\ 

finishes the proof. 

For the mixed moments we have by Lemma \2. 31 

E((TrU n ) k (TF(j;r) =0 (k^m), 

and we have proven the convergence of all moments. The only thing is left to conclude 
the convergence in distribution is to show that the moments determine uniquely the 
limiting distribution ( VIII. 6 in [I]). By the Carleman criterion (in VII. 3 of for a 
real valued random variable the moments jk determine the distribution uniquely if 



oo. 

ken 



Although we have complex random variables, the distribution of the argument is uni- 
form, and we can consider them as real valued random variables. The Stirling formula 
tells us 



k&i k>M \ v 7 / k>M 

for a large M G N, since V2k~K < 2 fc , if k > 2. □ 
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4 Asymptotic behaviour of the traces of higher pow- 
ers 



The aim of this section is to study the trace of higher powers of a Haar unitary. This 
was done also by Diaconis and Shashahani in Here we use elementary methods. 

Theorem 4.1 Let Z be standard complex normal distributed random variable, then for 
the sequence of U n n x n Haar unitary random matrices Tr U l n — > ylZ in distribution. 

Proof. We use the method of moments again. Lemma 12.31 implies that we only have 
to take into consideration E ^(Tr U l n ) k ^Tr U^j J , for all k E Z + . 



E {(TrU l n ) k (TrUi] 



3l,—,3l 



u, 



l I(fc-l) + l l i(k-l)+2 



u, 



*fci l !(fe-l) + l 



x Z7j 13 - 2 . . . U jtf-JJ j i+1 ji +2 ■ ■ ■ V j2ih+i • • • Uji( k -i)+di(k-i)+2 • ■ ■ "ifciii(fc-i)+i / ' 

where the indices h, . . . , ikhji, ■ ■ ■ ,jki run from 1 to n, and by Lemma 12731 if the sets 
{ii, . . . , iki} and {ji, . . . ,jki} are different, then the expectation of the product is zero. 
It follows from the the Cauchy and Holder inequalities, and (jHJ), that 



E yUi^ . . . Ui kl i l(k _ 1)+1 U j 1 j 2 . . . U j kl j l(k _ 1)+l 



< E 



U^iz . . . Ui kt i l(k _ 1)+1 U j 1 j 2 . . . Uj kl j [{k _ 1)+1 



< \lE(\U ili2 



I 2 \U- ■ 

I • • • \^ ^klH(k-l)+l 



2 \U ■ ■ I 2 \u 



3kl3l(k- 



, )+l ?)<0(n- kl ). 



Again the number of the set of indices, where there exist at least two equal indices is at 
most 0(n kl ~ r ), so the sum of the corresponding expectations tends to zero as n — > oo. 
Suppose that all the indices are different. There exist 



(n-kl)\ 



(kl)\ = 0(n M ) of these 

kinds of index sets, and now we will prove, that most of the corresponding products 
have order of magnitude less than n~ kl ^ 1 . Consider for any < r < kl 



NJt(r) := E (\U 12 \ 



1 23 1 2 • • • \U r l\ 2 U r+ i t r + 2 ■ ■ ■ Ukl-l,klUkl,r+\U r+2,r+l • • • U r +l,kl) ■ 



Note that NJ?(kl) = N%(kl - 1) = M%, and if {i u . ..i u ] = {j u . . . ,j H }, and all the 
indices are different, then the corresponding term equals to N£(r) for some < r < kl. 
Using the orthogonality of the rows for < r < kl — 2 



M5> 



12 



\u 23 \ 2 ...\u rl \ 2 u, 



r+l,r+2 



• ■ ■ Ukl-ljU, 



,j<Jkl,r+lU r+2,r+l 



u 



r+l,j 



0. (12) 
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If j > kl, then the permutation invariance implies, that 

E ( | 12 | 2 1 23 1 2 • • • \U r l\ 2 U r+ i <r+ 2 ■ ■ ■ Ukl-l,jUkl,r+\U r +2,r+l ■ ■ ■ Ur+l,j) = N£(r), 

so we can write from (|12p 



kl 

' I r r_ _ 1 2 i r r 1 2 \tt a 2 tt . _ tl. , .tt.. .77 ._ 77 



{n-kl)N%(r) = -E ( ^ I^H^sT • • • |f/rl| 2 f/r+l,r+2 • • • U kl -ljU kljr+1 U r+2 , r+ l . . . U r+1 



On the right side there is a sum of kl numbers which are less than 0(n kl ) because of 
(HU, so this equation holds only if N£(r) < 0(n~ H ~ 1 ). 



We have to compute the sum of the expectations 

E (jUi^l ...If/j^J ■ ■ ■ |f^i (fc _ 1)i+1 i (fc _ 1)i+2 1 • ■ • \Ui kl i (k _ 1)l+1 1 j = M kl . 

Now we count the number of these summands, so first we fix the set of se- 
quences of length I Ii )k = {(i( u -i)i+i, ■ ■ ■ iKi), 1 < u < k}, and we try to find the set 
Ji,k = {{j(u-i)i+i, ■ ■ - Jul), 1 <u < k}, which gives Mfo. If the product contains U irir+1 , 
then it has to contain Ui r i r+1 , so if i r and i r+ i are in the same sequence of Ii sk , then 
j s = i r and j t = i r+ \ have to be in the same sequence of Ji tk , and t = s + 1 modulo I. 
This means, that for all 1 < u < k there exists a sequence . . . , i v {) G I k ,i and a 

cyclic permutation n of the numbers {(v — 1)1+1, . . . , vl} such that • • • ,j u i) — 

(i7r((«-i)j+i) , • • • )«7r(«o)- We conclude, that for each 7^ there are Ji jk , since we can 
permute the sets of I^ k in k\ ways, and in all sets there are I cyclic permutations. 
Clearly there are ^' kl ^ sets Ii jk , so 

lim E ( (TrUtf (t^Y) = lim , n \ kll ^ ~ £ ^ ~ k ^ = k \l\ 
and as in the proof of Theorem 13. II this is the kth moment of (y/lZ)(y/lZ). □ 



5 Independence of the trace of the powers 

In this section we prove that the limits of the trace of different powers are independent. 
The method of computation is the same as in the previous sections. 

Theorem 5.1 Let U n be a sequence of Haar unitary random matrices as above. Then 
Tr U n , Tr U%, . . . , Tr U n are asymptotically independent. 

Proof. We will show, that the joint moments of Tr U n , Tr U%, ■ ■ ■ Tr U n converge to the 
joint moments of Z\, \J2Z2, • • • , \flZi, where Z\, Z2, ■ ■ ■ Z\ are independent standard 
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complex normal random variables. The latter joint moments are 




= U^ E { z " z ^) = Tl 6 ^ Hat 

1=1 i=l 

^From Lemma ESI if Yl\=i Yl\=i then the moment 

E[fl i^KY* (W) J =0- 



Now ^icii = 'Yl/ibi. Again if the indices in the first and the second part are not the 
same, then the expectation is zero according to Lemma \2. HI The order of magnitude of 
each summand is at most O as above, so if not all the indices are different, 

then the sum of these expectations tends to zero, as n — > oo. The same way as in the 
proof of the previous theorem, those summands where there is a U irir+1 U iria1 i r+ \ ^ i s 
in the product are small. So now we have to sum the expectations Mj^ ia .. If we fix the 
set of first indices /, then again the sequences of the appropriate J, have to be cyclic 
permutations of the sequences of I. So the number of the sequences of length i in / is 
the same as in J, which means a, = 6j for all 1 < i < I. The number of the / sets is 
t — w . so we have arrived to 

(n-^ ia,i)\ ' 



n ^°° \i=i 



lira =— HS^aA = \[8 aifii aS ai ■ 

i^oo [n — > id; ! - LJ - n! - LJ - 



□ 



6 Truncation 

Let U be an n x n Haar distributed unitary matrix. By truncating n — m bottom rows 
and n — m last columns, we get an m x m matrix U\ n , m ]- In this section we study the 
limit of U[ n>m ] when n — > oo and m is fixed. Our method is based on the explicite form 
of the joint eigenvalue density. 

The truncated matrix is not unitary but it is a contraction. Hence the eigenvalues 
Z\,Z2, ■ ■ ■ , z m G D m , where D = {z G C : \z\ < 1} is the unit disc. According to 19, 
the joint probability density of the eigenvalues is 

rn 

i<j i=l 
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on D m . 

Since the normalizing constant C\ n , m \ was not given in PJ, we first compute it by 
integration. To do this, we write Q = Tie 1<Pi and dQ = r, dri d(p{. Then 



n i^-^nc 1 



|2\n— m— 1 



l<i<j<m 



8=1 



[0.1] m - / [0^r] m i<2< r 



8=1 



1=1 



Next we integrate with respect to dip = d<pi dtpi . . . d(p m by transformation into com- 
plex contour integral what we evaluate by means of the residue theorem. 

[°' 2 ^ n l<i<j<m 



n 

l<i<j<m 



i=i 

n - r J%)( r ^r i - r i z j i ) n dz 

' n Ki<i<m i=l 



„ m 

/ n^ idet 



1 

T\Z X 



1 

^2 



m-1 m-1 ra-1 m-1 



2 ^2 



x det 



riZi 



^2 



m-1 m-1 -(m-1) 

1 1 z 1 '2 z 2 



r m,Z m 



r m—l m-1 



T rn z m 



„m-l~-(m-l) 

I m Z m 



X 



„ m m m 

= nr / IR 1 E (- 1 )* 00 II^ " 1 E Ifr* 1 )* 0-1 

^ i=l 7rG5 m i=l pG5 m 8=1 

We have to find the coefficient of YllLi 3 this gives that only p = it contribute and 
the integral is 

(2vr) m £ n^) 2(PW_1) - 

pe5 m i=l 

So we have 

„ m m m 

q-i = (2-) m / £ nN 2(p(iM) - ^~ m ~ 1 n ^ ^ 

■ / [°. 1 l m pes™ <=i «=i i=l 

m »i 

= (27r) m m!n / ^(l-r?) 



8=1 



2\n— m— 1 



dr. 
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and the rest is done by integration by parts: 



/ 

Jo 



„2fc+l, 



1 — r 



2\n— m— 1 



dr 



n — m 



2k-l, 



1 — r 



2\n— m 



dr 



k\ 



n — m) . . . (n — m + k — 1) J 

-i 



r(l — r ) 



2\n— m+fc— 1 



n — m + k — 1 
k 



1 



2(n — m + k) 



Therefore 



ir m m 



m-l / , i \ -1 

j-j- / // — /// /. — I 



k=0 



k 



n — m + k 



Now we consider y/n/mU\ njm ]. Its joint probability dnsity of the eigenvalues is simply 
derived from the above density of C/[ n>m ] by the transformation 



(G> •••■> Cm) 



m m 

Sli ■••it/ ^m I ? 

n V n 



and it is given as 



n — m + k — 1 
k 



i<j v 1 ; i v 

1 m— ] 

= — TT 

7r m m! 11 

fc=0 - 

xnic-c/nfi-^) 

i<7 i=l ^ 

1 m— 1 

= — TT 

7r m m! 11 

fc=0 



2 \ n—m—1 



(n — m + fc) ^ — j 

2 \ n—m—1 



m(m+l)/2 



n *+ 1 (l + o(l)) /m\">Ki)/ 



II III 



2 \ n—m—1 



rn m(m+l)/2 m / ttiI^I 2 \ 

^( 1+0(1 ),nio-c/n0-^) 



i=i 

2 \ n—m—1 



The limit of the above as n — > oo is 

m m(m+l)/2 



exp -m^|0| 2 J II IC* - C?r» 



i=i 



which is exactly the joint eigenvalue density of the standard m x m non-selfadjoint 
Gaussian matrix. This implies the following theorem. 
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Theorem 6.1 The normalized truncated matrix 

U[n,m] 

converge in distribution to the standard m x m non- self adjoint Gaussian matrix. 
Details of this convergence will be subject of a forthcoming publication. 
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